Symmetry causes a huge conductance peak in double quantum dots. 
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We predict a huge interference effect contributing to the conductance through large ultra-clean 
quantum dots of chaotic shape. When a double-dot structure is made such that the dots are the 
mirror-image of each other, constructive interference can make a tunnel barrier located on the sym- 
metry axis effectively transparent. We show (via theoretical analysis and numerical simulation) that 
this effect can be orders of magnitude larger than the well-known universal conductance fluctuations 
and weak-localization (both less than a conductance quantum). A small magnetic field destroys the 
effect, massively reducing the double-dot conductance; thus a magnetic field detector is obtained, 
with a similar sensitivity to a SQUID, but requiring no superconductors. 
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In the 1990s, interference effects (universal conduc- 
tance fluctuations and weak-localization) were observed 
for electrons flowing through clean quantum dots [l|, Q ■ 
The chaotic shape of such dots makes these effects analo- 
gous to speckle-patterns in optics rather than to the reg- 
ular interference patterns observed with Young's slits or 
Fabry-Perot etalons. While such interference phenom- 
ena are beautiful, they have only a small effect on the 
properties of quantum dots coupled to multi-mode leads. 
Here we provide a theoretical analysis and numerical sim- 
ulations showing that a much larger interference effect 
occurs in systems which are mirror-symmetric but oth- 
erwise chaotic 0, i, i, @, see Fig. □ We show that the 
mirror symmetry induces interference that greatly en- 
hances tunneling through a barrier located on the sym- 
metry axis; it can make the barrier become effectively 
transparent. Thus an open double-dot system with an 
almost opaque tunnel barrier between the two dots will 
exhibit a huge peak in conductance when the two dots 
are the mirror image of each other, see Fig. [51 This effect 
could be used to detect anything which breaks the mirror 
symmetry. For example, current 2D electron gas (2DEG) 
technology 0] could be used to construct a device whose 
resistance changes by a factor of ten, when an applied 
magnetic flux changes from zero to a fraction of a flux 
quantum in the double dot. This is a sensitivity similar 
to that of a SQUID, but it is achieved without supercon- 
ductivity, making it easy to integrate with other 2DEG 
circuitry. 

Origin of the conductance peak. The origin of the 
effect can be intuitively understood by looking at Fig. [T] 
Assume that electrons only follow the two paths shown 
(instead of an inflnite number of different paths). Path 1 
does not tunnel the flrst time it hits the barrier, but does 
tunnel the second time it hits it. Path 2 tunnels the first 
time it hits the barrier, but not the second time. Quan- 
tum mechanics gives the probability to go from the left 
lead to the right lead as \r{e)t{0')e'^^/'^+t{e)r{e')&'^^/'''\^, 
where the scattering matrix of the tunnel barrier has am- 
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FIG. 1: A mirror-symmetric double dot, where the classical 
dynamics is highly chaotic. We call it a "butterfly double 
dot" to emphasize the left-right symmetry. Every classical 
path from the left lead to the right lead (solid line) which hits 
the barrier more than once, is part of a family of paths which 
are related to it by the mirror symmetry (dashed line). 

plitudes r{9) and t{9) for reflection and transmission at 
angle 6. If there is no correlation between the classical 
actions of the two paths {Si and S2), then the cross-term 
cancels upon averaging over energy, leaving the probabil- 
ity as \r{9)t{e')\'^ + \t{e)r{e')\'^. In contrast, if there is a 
perfect mirror symmetry, then S2 = Si, and the probabil- 
ity is \r{9)t{9') + t{9)r{9')\'^ , which is significantly greater 
than |r(0)i(V)|2 + \t{9)r{9')\^ . Indeed, if we could drop 
the 6'-dependence of r and t, the probability would be 
doubled by the constructive interference induced by the 
mirror symmetry. A path that hits the barrier (n -|- 1) 
times has 2" partners with the same classical action (each 
path segment that begins and ends on the barrier can 
be reflected with respect to the barrier axis). However 
the conductance is not thereby enhanced by 2", because 
(due to the nature of the barrier scattering matrix) there 
is also destructive interference when one path tunnels 
(4j — 2) times more than another (for integer j). 

The effect looks superficially like resonant tunneling. 
However, that only occurs when dots are weakly coupled 
to the leads, so that each dot has a peak for each level of 
the closed dot and the current flow is enhanced when two 
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FIG. 2: Average conductance as (a) a function of applied B- 
field (with the barrier on the symmetry axis), and as (b) a 
function of the barrier position (for zero B-field). The latter 
mimics the effect of gates that reduce the size of one dot 
relative to the other. The data points come from simulations 
performed for the structures shown in the insets. The curve 
comes from the semiclassical theory; in (b) there is no fitting 
parameter, while in (a) an unknown parameter (or order one) 
is adjusted to fit the data. The conductance of the tunnel 
barrier alone is Gtb. 



peaks are aligned. Instead in our case each dot is well 
coupled to a lead (with iV ^ 1 modes) , so the density of 
states in each dot is featureless (the broadenning of each 
level is about iV times the level-spacing). Furthermore, 
resonant tunneling occurs at discrete energies, while our 
effect is largely energy independent. Another similar ef- 
fect, called "reflectionless tunneling", occurs when elec- 
trons are retro-reflected as holes, due to Andreev reflec- 
tion from a superconductor [1, However, this retro- 
reflection transforms the classical dynamics in the dot 
from chaotic to integrable [loj , and large interference ef- 
fects in integrable systems are not uncommon (consider a 
Fabry-Perot etalon). Here, the mirror symmetry induces 
a large interference effect without any retro-reflection and 
without a change in the nature of the classical dynamics 
(chaotic motion remains chaotic). 

Semiclassical theory. Our analysis uses the semi- 
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FIG. 3: Plot of the ratio (Gsym)/(Gasym>, given by Eqs. (f5l6ll . 
The ratio grows as Ttb —> for all P (although (Gsym.asym) 
shrink). For given Ttb, the ratio is maximal at P — (1 — 

2r;/')/(i-4Ttb). 



classical theory of transport through clean chaotic quan- 
tum dots [ll| . The conductance through a system whose 
dimensions are much greater than a Fermi wavelength 
can be written as a double sum over classical paths, 7 
and 7', which both start at a point yo on the cross-section 
of the left lead and end at y on the right lead: 

G = (27rfi)-iGo^A^A;,exp[i(5-y-5y)/n], (1) 

7,7' 

where Gq ~ 2e'^/h is the quantum of conductance, and 
S-y is the classical action of path 7. A tunnel barrier with 
left-right symmetry must have the scattering matrix 
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where r(9) and t(9) are reflection and transmission am- 
plitudes for a plane wave at angle of incidence 6. Keep- 
ing only the upper sign in 5tb(^) fl3| , the amplitudes in 
Eq. H]) are 
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where path 7 starts with a momentum across the left 
lead, pyg, and a total momentum given by the Fermi mo- 
mentum, pp. This path reflects off the barrier mn^j) 
times (with the fcth reflection at angle 9B.k) and trans- 
mits mT(7) times (at angles 9Tk) before hitting the right 
lead at y. The factor {dpyg/dy)j is the stability of the 
path that would exist if the barrier were absent for each 
transmission and impenetrable for each reflection. For 
most pairs with 7 7^ 7', the exponent in Eq. ([1]) varies 
fast with energy, so that averaging over energy removes 
such pairs from the double sum. We keep only the main 
contributions surviving such averaging: those where 7' 
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can be constructed from 7 by means of the reflection with 
respect to the barrier axis (symmetry axis) of any path 
segment that begins and ends on the barrier, for which 
*S''y' ^ at all energies (the paths thereby have the same 
stability (dpyg/dy)^). Dropping weak-localization effects 
the average conductance reads 



dy 



"(7) 



(4) 



where the product is ordered, and n{'-f) is the number 
of times the path 7 hits the barrier. The four-by-four 
matrix § = iStb ® S^^ gives the scattering of the two 
paths at the barrier. Thus S^- gives the weight to go 
from state j to state i, where we define state 1 as both 
paths in the left dot; state 2 as path 7 in the left dot and 
path 7' in the right dot; state 3 as path 7 in the right dot 
and path 7' in the left dot; and state 4 as both paths in 
the right dot. The matrices C^^m are diagonal with the 
following non-zero elements: [C-,..m]ii = [C-y,m]44 = 1 
and [C-y_rri]22 = [C^^mJas = exp[i(55-y_m/fi-]- The action 
difference SS^^m is that between path 7 in the left dot and 
its mirror image in the right dot between the (m — l)th 
and mth collision with the barrier. For perfect symmetry 
C„i = I and then the product equals [§"]4i. 

We assume that the classical dynamics is sufficiently 
mixing that paths uniformly explore the dot between 
subsequent collisions with the barrier (or leads). Defin- 
ing SSq/H as the phase difference acquired in one time of 
flight across the dot, we have C^^m — cxp[— Ft^ where 
F is a complex number, with Im[F] ~ (SSq) /{toH) and 
Rc[F] ~ var[(5S'o]/(To?i^). The probability that a path 
survives in the dot for a time t without hitting either the 
barrier or the lead is e"*/'^^ . Using this, we replace C^.m 
by its time-average C — (C-y,,,,.); its only non-zero ele- 
ments are Cn = C44 = 1 and C22 = C^g ^ [I + Fr^]-^ 
Thus the product in Eq. ^ reduces to (C§)". The sum is 
over all 7s that hit the barrier n times, and is independent 
of yo, y- To proceed, we define S = C^/^SC^/^; it is simple 
to show that [(CS)"]^^ = [§"]4i for aU n. Then, defining 
P Wtb/(W^tb + W) as the probabihty for a path to hit 
the VFtb-wide barrier before escaping into the VF-wide 
lead, we find that (G) = GoiV(l - ^') E^^i ^" [^"] 41' 
where N = p-pW/{T:h) is the number of modes in a lead. 
Upon finding the matrix, U, which diagonalizes S, one 
can easily evaluate the geometric series in n. 

This analysis gives the following average conductance 
of the symmetric double dot (F = 0), 

(Gsym) = GoiVP(l + P)Ttb/[(l-P)'+4PTtb], (5) 

where Ttb is the tunneling probability, \t{9)\'^, averaged 
over all 9. For Ttb < (1 - P)/2 (i.e. for Gtb, the conduc- 
tance of a barrier with transmission Ttb in a waveguide 
of width Wtb, less than P times the conductance of the 
scries of the two constrictions), one finds that (Ggym) is 



greater (often much greater) than the tunnel barrier con- 
ductance, Gtb- Thus symmetrically placing constrictions 
on either side of the barrier can strongly enhance its con- 
ductance (this is a stark example of the fact that quan- 
tum conductances in series are not additive). In contrast, 
for the asymmetric double dot (large F) we have 



(Ga.ym) - G„NPT,^/[l-P+2PTt 



tbj 
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which is always less than Gtb- The ratio (Gsym)/ (Gasym) 
is plotted in Fig. [3] For any finite Ttb, the ratio is max- 
imal at P = (1 - 2Tj^/^)/(l - 4rtb). This choice of P 

gives (Gsym) = Go^/4 and (for small Ttb) (Gasym) ^ 
1/2 

T^^^ GqN/2. Thus the conductance ratio can be arbi- 
trarily large for a highly opaque tunnel barrier. 

Peak shape with symmetry-breaking. The effect 
of the mirror symmetry is suppressed by a perpendicular 
magnetic field, i?, or by moving the boundary of one 
dot by a distance 5L. It is also suppressed by disorder 
(defined by a mean free flight time between subsequent 
scatters from disorder, Tmf ) or decoherence (defined by a 
decoherence time, r^). The suppression can be quantified 
in terms of the following parameters: 
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where e is the electronic charge, A is the area of one dot, 
and To is the time to cross the dot. In Fb, the constant r] 
is of order one, but is hard to estimate [l^l- For Fboundary, 
we have {5L) ^ x£, and va.r[dL] ~ a;^(^ — ^^), if a fraction 
^ of the left dot is deformed outwards by a distance x. 
For multiple asymmetries, the total F is the sum of the 
individual Fs given above. For real F, 



(G(F)) = (Gasym) + 



{Gsym ) (Gasym ) 

l + F(P,Ttb) xFr/3' 



(10) 



where T(P,Ttb) = (Gsym)/[(Gasym)(l + P)], and ~ 
■7tLto/{W + Wtb) is the typical time a path spends in 
one dot before hitting a lead or the barrier. For the 
large conductance ratio (see below Eq. F{P, Ttb)TQ 
is about half the dwell time in the double-dot, td ^ (1 — 
P)~^T^. Thus the conductance is a Lorentzian function 
of the P-field, with similar width to the weak-localization 
dip in the same system with no barrier IJ]. This makes 
the system an extremely sensitive detector of magnetic 
fields and deformations of the confining potential (for 
example due to the movement of charge near the double 
dot). Intriguingly, the peak remains when the leads are at 
different positions on the two dots; it is simply suppressed 



with an asymmetry parameter F 



lead 



(l-P)/7 



For complex F, as in Fig. Eljb), we have no analytic 
result for (G(F)), but we can get it by numerically di- 
agonalizing the 4-by-4 matrix, S. In Fig. [DJb), the data 
and the theory curve drop below (Gasym) — 0.23Go. We 
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will show elsewhere that this is due to destructive inter- 
ference. The conductance rises back up to (Gasym) when 
the barrier is moved a distance of order a wavelength. 

Proposal for experimental observation. Con- 
sider making such a double-dot in an ultra-clean two- 
dimensional electron gas (2DEG) at the lowest achiev- 
able ternperatures 0. A finger gate could define the 
barrier [l6|, with split gates controlling the lead widths. 
To maximize the effect for a 2DEG with a mean free 
path 0] of order 500 fim, each dot (see Fig. 2) can have 
size L = 4 /im (circumference ~ 3.6L ~ 15 ^m) with 
12 mode leads {W = 310 nm ~ 6Af). A barrier with 
Ttb = 1.48 X 10^3 and width Wth = L gives P = 0.93 and 
r/) ^ 3.5ro. In this case, (Gsym) ^ 14(Gasym) ^ 3.2Go 
(resistance i?sym ^ 5kQ). The crossover from (Gsym) 
to (Gasym) happens for T ~ 0.14/rD 0.04/to. At low 
temperatures {t^ > Tmf), disorder will suppress the peak 
to about 83% of (Gsym), since i^(P, Ttb)r„fT/) - 0.2. 
Thus the double-dot conductance will drop by an order 
of magnitude if 10% of the boundary of one dot is moved 
by Af/2, or if a B-field is applied such that a fifth of a 
flux-quantum threads each dot. The latter is a B-ficld 
sensitivity similar to that of a SQUID. 

The main experimental challenge will be to define dots 
that are mirror-symmetric on a scale significantly less 
than Af ^ 50 nm. We suggest that each dot should be 
defined by means of multiple gates (made as symmetric 
as possible); their voltages can then be tuned to max- 
imise the symmetry. We propose the following protocol 
for this maximization. Starting with very wide leads, 
in such a way that P is far from unity and the conduc- 
tance peak is very broad, one scans the dot-defining gate 
voltages over a broad range to reveal the approximate 
symmetry point (maximal conductance). One then nar- 
rows the leads (increasing P), making the conductance 
peak higher and narrower, and adjusts the dot-defining 
gate voltages to again maximize the conductance. Re- 
peatedly doing this should give the symmetry point with 
increasing accuracy, until one reaches the limit imposed 
by inherent asymmetries (disorder, etc). 

Numerical simulations. For the above maximisa- 
tion we took Wth = L and only 12 lead modes. This 
calls into question two assumptions in the theory. Firstly, 
we can no longer assume that paths in the dot will be 
well randomized between collisions with the barrier, since 
Tq ~ 3.6ro. Secondly, we may not be able to neglect 
other interference effects (weak-localization, etc), since 
(G) is at most a few Gq. Thus to verify that the effect 
is as expected in such a parameter regime, we numer- 
ically simulated a stadium billiard containing a barrier 
with Ttb = 1-48 X 10~^, see Fig . [2] We use the recursive 
Green's function technique [17| working in real space for 
the direction of current propagation (cut into multiple 
slices) and in mode space for the transverse direction. 
Magnetic fields are in a Landau gauge where the vector 
potential is oriented in the transverse direction [isj . The 



number of longitudinal slices and transverse modes were 
increased until the results converged. The data shown 
here are for 836 longitudinal slices (200 of which are in 
the outer leads) and 200 transverse modes. We mimic 
thermal smearing, at a temperature of 23 mK, by aver- 
aging over 44 energies uniformly distributed over an in- 
terval of 0.02 meV around the Fermi energy of 9.02 meV. 
We use the effective mass in GaAs of 0.067 mg. The simu- 
lation (data points in Fig. [2]) clearly shows that the effect 
exists in this regime. Indeed, despite the assumptions in 
its derivation, the theory (solid curve) agrees surprisingly 
well with the numerical data. 

Concluding comment. The conductance peak is not 
destroyed by bias voltages or temperatures greater than 
H/to, because the mirror symmetry is present at all en- 
ergies and not just at the chemical potential (unlike the 
electron-hole symmetry for reflectionless tunneling into 
a superconductor). Large biases or temperatures should 
still be avoided, as they increase the decoherence. We 
thank M. Houzct and P. Brouwer for discussions. 
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APPENDIX (ONLY ON ARXIV VERSION) 

Comment on the conductance ratio. It is instruc- 
tive to consider particular limits of Eqs. (j5|6p . At the 
symmetry point we observe that the more times a path 
returns to the barrier, the more transparent the inter- 
ference makes the barrier. Thus if the path takes an 
infinite time to escape the double dot {P ~ 1), then the 
barrier becomes completely transparent. However this 
does not generate a large conductance peak, since for 
P — > 1 the probability to go from the left lead to the 
right lead is a half for any finite barrier transparency, 
thus (Gsym)/(G asym) — * 1- This is the reason for the 
maximum conductance peak occurring when P is slightly 
less than one (as visible in Fig. [3]). The opposite limit 
is P ^ 0, then both (Gsym) and (Gasym) reduce to the 
conductance of the barrier alone. Since no path hits the 
barrier more than once, there can be no interference in- 
duced enhancement of tunneling. 

Comment on fitting i?-field dependence. To 
make the theory quantitative (for comparison with the 
numerical simulations) we assumed that the area en- 
closed by each straight path segment from one point 
on the boundary of the left dot to another is uncorre- 
lated with the next. We define this directed area, a, as 
that of the triangle made by the two ends of the path 
segment and the mid-point of the barrier. Assuming 
{eBa/h) ^ 1 we have Tb = 2K{eBA/hy /tq, where the 
system-specific parameter k = var(a/^). If a were uni- 
formly distributed over the range from — A/ 2 to A/ 2, we 
would get K = 1/12. In contrast, if the distribution were 
strongly peaked at —A/ 2 and A/ 2, we would get k as big 
as 1/2. We believe only a ray-tracing simulation of the 
cavity would yield an accurate value for k, thus for the 
theory curve in Fig. [2^ we treat At as a fitting parameter. 
The Lorentzian width of 40 ^T corresponds to k = 0.17; 
which is within the range estimated above. 

Comment on fitting barrier-position depen- 
dence. In the case shown in Fig. Oa, the mirror sym- 
metry is only broken at the point where a path segment 
begins or end at the tunnel barrier. Thus a path segment 
only acquires a phase-difference from its mirror image at 
the places where it touches the barrier. Paths acquire 
more phase in the left dot than the right dot (since the 
barrier is moved to the right). This phase difference has a 
very different form from that induced by a i?-field (where 
the phase difference grows with the time a path spends 
in one of the dots). Taking this into account, we get the 
solid curve in Fig. Wijo) without any fitting parameters. 



